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The propagation of light in a resonant atomic vapor can a priori be thought of as a multiple 
scattering process, in which each scattering event redistributes both the direction and the frequency 
of the photons. Particularly, the frequency redistribution may result in Levy flights of photons, 
directly affecting the transport properties of light in a resonant atomic vapor and turning this 
propagation into a superdifusion process. Here, we report on a Monte-Carlo simulation developed 
to study the evolution of the spectrum of the light in a resonant thermal vapor. We observe 
the gradual change of the spectrum and its convergence towards a regime of Complete Frequency 
Redistribution as the number of scattering events increases. We also analyse the probability density 
function of the step length of photons between emissions and reabsorptions in the vapor, which 
governs the statistics of the light diffusion. We observe two different regime in the light transport: 
superdiffusive when the vapor is excited near the line center and normal diffusion for excitation far 
from the line center. The regime of Complete Frequency Redistribution is not reached for excitation 
far from resonance even after many absorption/reemission cycles due to correlations between emitted 


and absorbed frequencies. 


PACS numbers: 42.25.Dd,05.40.Fb,42.68.Ay,32.70.Jz 
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I. INTRODUCTION 


Light scattering in atomic vapors has been studied for a long time, for resonant as well as for non-resonant 
radiation l|, When light is resonant with atomic transitions, many absorption-reemission cycles may 
occur before a photon leaves the vapor volume, a phenomenon known as imprisonment of resonant radiation 
3| or radiation trapping This multiple scattering process alters the spatial distribution of excitation in 
a resonant vapor, as well as the frequency spectrum of the light, because photons usually suffer a frequency 
shift at each scattering event. The importance of the frequency redistribution for light diffused in resonant 
vapors was already acknowledged in the 1930s. Radiation trapping and frequency redistribution are long- 
studied topics in astrophysics, for instance in the analysis of the radiation emitted by nebulae 5l-l8|. 


The broadening of spectral lines favors the emergence of long steps for photons in a resonant vapor 


such a way that the light escapes the vapor volume faster than expected in a diffusive description 
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Occasional very long steps of photons result in the impossibility of defining a mean square displacement. 
To overcome the failure of a diffusive model to describe light transport in resonant vapors, Holstein 3| 
developed a Boltzmann-like integro-differential equation taking into account the strong spectral variations 
of the absorption coefficient around resonance frequencies. Holstein equation is usually solved for Complete 
Frequency Redistribution (CFR), i.e. assuming that there is no correlation between the frequency of the 
absorbed (incident) and of the reemitted photons. In this case, the probability density function (PDF) 0(i/) 
of the incident radiation is equal to the PDF $(i^) of the vapor absorption l^ll3|. CFR usually occurs when 
atomic collisions are very frequent and destroy the correlations between incident and scattered photons [7 1. 


For low-density vapors, where interatomic collisions can be ignored, the frequency of the reemitted photon 
is exactly the same as the frequency of the absorbed photon, in the atom’s rest frame (ARF), i.e. there 
is no frequency redistribution in the ARF. However, in the referential of the laboratory, the Doppler effect 
due to the motion of the scattering atoms shifts the frequency of the photons, leading to Partial Frequency 
Redistribution (PFR) in the laboratory frame. One can define a redistribution function R{v'^v) that gives 
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the probability R(y'^v)dv'dv that a photon of frequency between u and u + dv \s absorbed and that the 


reemitted photon has a frequency between v' and v' + dv'. This redistribution function was calculated by 
Unno 


for an infinitely sharp but Doppler-broadened line. Latter, Unno 


14j and Hummer 


15l| calculated 


the redistribution function for a line profile with both natural and Doppler broadening. As Holstein equation 
is not easy to solve in PFR conditions, the study of light transport in resonant scattering media is usually 


carried out through Monte-Carlo (MC) simulations 


id. Such numerical simulations allow to follow both 


the spatial distribution of the light excitation and the frequency spectrum of the light after a given number 
of scattering events. MC simulations have been used to obtain information such as the number of scattering 
events that a photon undergoes before it escapes from a vapor cell with a particular geometry 
local thermodynamic equilibrium assumed by Holstein 17|; study frequency diffusion when a cloud of cold 


161; test the 


atoms is illuminated by a laser light detuned up to a few natural widths from resonance 


18 l. 


Besides its applications in astrophysics, radiation trapping is also a phenomenon of great interest in 
atomic physics laboratories. Radiation trapping increases the effective lifetime of the excited population 
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Radiation trapping diffuses light incoherently and contributes to ground state decoherence in experi¬ 


ments of Coherent Population Trapping (CPT) and Electromagnetically Induced Transparency (EIT) 


The diffusion of photons in a vapor modifies the spatial distribution of excited atoms 
role in the spatial distribution of saturation 221. 
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21j and thus plays a 


The strong dependence of the transport properties of light on the spectral characteristics of the atom- 
photon interaction has led to an increasing interest in studying the statistics of the random walk of photons 
in a resonant vapor. It has been theoretically predicted that tlm length of the steps between two scattering 
events does not follow the statistics of normal distributions 
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25l| . Anomalous, superdiffusion of light in 


atomic vapors occurs due to long steps taken by reemitted photons whose frequency lays far in the aisles 
of the curve of the absorption probability distribution. Superdiffusion is then characterized by a PDE of 
step length I that asymptotically follows a power law P{1) ~ l~^ with fc < 3. The step-length distribution 
depends on the spectrum of the incident radiation as well as on the absorption profile of the scattering 
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vapor 2^. All the standard atomic spectral lineshapes (Doppler, Lorentzian and Voigt) result in such a 
superdiffusive behavior, characterized by rare, large steps known as Levy flights. For Voigt incident and 
absorption spectra, for instance, the theoretical prediction yields 23j k = 3/2. While it is trivial to observe 


the exponentially decaying step-length distribution of (monochromatic) laser photons in a Voigt-broadened 
vapor, measuring the PDF of the photons step length in the CFR regime is not an easy task. To do so, one 


may, for instance, prepare an incident radiation with a Voigt-like profile, which was achieved in Ref. 


26| by 


submitting laser photons to a few frequency-redistributing scattering events in an auxiliary vapor cell prior 


to sending them to the measurement cell 


27|. 


Here we are interested in this ’’preparation” of a Voigt spectrum for the radiation for studies of photons 
steps statistics. We implement a MC routine based on first principles, whose aim is to analyse the PFR in 
the few initial scattering events. In particular, we want to determine how many scattering events are needed 
to generate an emission spectrum similar to the absorption one, i.e., a CFR. From this simulation we can also 
obtain the evolution of the PDF of the photons step length. The knowledge of the incident spectral profile 
and of the PDF of the step size is crucial in the experimental study of the statistics of photon diffusion in a 
vapor 26, 27|. 


II. MONTE-CARLO SIMULATIONS 

In this section we describe the Monte-Carlo simulation used to study the frequency redistribution of 
photons scattered in an atomic vapor. 

We consider a photon in a resonant vapor and departing from the origin of the coordinate system with 
a detuning 6i relative to the center of the atomic transition and with an isotropic random direction. This 
photon may be absorbed by atoms of the vapor and then reemitted many times. We call i the scattering 
event number. After the ith scattering event the photon is reemitted isotropically, with a detuning 6i related 
to the absorbed frequency and to the velocity of the emitting atom through Doppler effect. The photon 
travels a distance li in the vapor before being absorbed again and the process of reabsorption/reemission is 
repeated until a given number of scattering events is reached. The frequency of the photon and the velocity 
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of the emitting atom are recorded after each scattering event. Details of the Monte-Carlo simulation are 
given below. 


The MC simulation is developed for two-level atoms with natural broadening of their excited level and 
with Doppler broadening. For such a system, the absorption profile is a Voigt profile with Voigt parameter 
a = where F is the natural line width, F£> = 2u/A is the Doppler width at 1/e of its maximum. 


u = yj2kBT/m is the half width of the Maxwell-Boltzmann velocity distribution at temperature T, A is the 
wavelength of the atomic transition, m is the mass of the atom and ks is the Boltzmann constant. We use 
the value a = 0.01 corresponding to the D 2 transition of ®®Rb at room temperature (T = 294 K). At such a 
temperature, the pressure broadening is negligible (~ 1 kHz 2^) and will not be taken into account. 


A. Photon step length 


The distance h traveled by a photon between the ith emission and the {i + l)th absorption in the vapor is 
drawn from the Beer-Lambert probability distribution function, P{l)dl = a{Si) exp (—a{6i)l)dl, which gives 
the probability that a photon of frequency 5i be absorbed after having traveled a distance between I and 
l + dl see Refs.[l7|,l25j. a{Si) is the absorption coefficient at detuning 5i as measured in the laboratory frame. 
It is proportional to the Voigt profile and to the atomic vapor density. 


B. Velocity of the absorbing atoms 

An important variable in the determination of the frequency redistribution by radiation trapping is the 
velocity of the atoms that absorb and reemit the photons. The speed of the atoms in a vapor follows a 
Maxwell-Boltzmann (MB) distribution. However, for a given frequency profile of the incident radiation, 
some atoms will absorb the photons more favorably, depending on their velocity component in the photon 
direction. The probability R 5 i(V||) that an atom with a parallel component of velocity V|| (positive sign in 
the direction of the photon) absorb a photon with detuning Si is given by the product of the Doppler-shifted 
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Lorentzian atomic lineshape, centered at ASi, and the Maxwell-Boltzmann distribution of Vj|: 


where SA,i = — (V||/A) is the detuning of the incoming photon in the atomic rest frame and (Vj|/A) is 

the Doppler shift. The probability density function of V|| of the absorbing atoms is shown in Figured] for 
different detunings of the incoming photons. 


The normal component of the atom velocity plays no role in the absorption process and is therefore 
independently drawn from a Maxwell-Boltzmann distribution at the temperature of the vapor. 

The PDF of the parallel velocity component is composed of two physically distinct contributions (Equation 
d]), with the wide MB distribution modulating the amplitude of the Doppler-shifted atomic Lorentzian 
lineshape. Due to this convolution of effects, the PDF of V|| (Figure 1) evidences two qualitatively different 
situations: i) for an incident frequency close to the atomic resonance (small Si), the Maxwell-Boltzmann 
probability of finding an atom with parallel velocity Vj| = ASi is high. The PDF of Vj| is a Lorentzian 
peak around V|| = ASi. This situation corresponds to an absorption at line center in the atomic rest frame 
(SA,i = 0, with a width AV|| AF, see Figures dla) and Figures dlb)(c)) for Sj = —^Fd and Sj = — Fd, 
respectively. Note that for Sj = —Fd (Figure die)) the photons’ detuning is larger than the Doppler width 
and the MB probability of finding an atom with parallel velocity V|| = AS/ is one order of magnitude 
smaller than at resonance. The probability distribution of the parallel velocity component is still essentially 
a Lorentzian peak around Vj| = AS/ (P(0) < P{ASi)) but atoms absorbing in the wings of this Lorentzian 
peak begin to give a noticeable contribution (broad peak at right of the narrow Lorentzian one), ii) For an 
incident frequency far from resonance, the probability of finding an atom with parallel velocity V|| = ASi is 
very small, because of the very fast decay (e ii) of the Maxwell-Boltzmann distribution. The photon is 
preferentially absorbed by an atom moving relatively slowly in the direction parallel to the incoming photon, 
i.e. inside the MB width but in the wings of the Lorentzian line-shape in the atomic rest frame (see Figures 
dla) and die)). This situation occurs when , that is, when the probability P(0) that 

the absorbing atom has V|| = 0 is higher than the probability P{ASi) that its parallel component of velocity 
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FIG. 1: (Color online) Probability density function of the velocity component of absorbing atoms parallel to the direction of 
incident photons, for different excitation detunings Sj. The values of Vjj are normalized to the half width of the 
Maxwell-Boltzmann distribution at temperature of 294 K. (a) In log scale, (b-e) in linear scale (zoom), (b) for Sj = 
and (c) for Sj = —Td: Photons with a small detuning are preferentially absorbed by atoms whose parallel velocity component 
Doppler compensates for the detuning, (d) for 5j = — I.ZPd and (e) for 5j = — 2r£): very few atoms absorb the far-detuned 
photons at the line center in the atomic rest frame, see small narrow peaks. Most atoms absorb instead in the far-reaching 
wings of the Lorentz lineshape, where favorable MB probability compensates for the very weak absorption probability (this 

one decaying as in the wings). 


is V|| = XSi. 

We will call the detuning corresponding to the limit between those two behaviors for which P(0) = 
P{Xdi). For the parameter a = 0.01 considered here, the limit between the two situations is around 
5i = ±6 l = il-TT^i (see FiguresIHa) and[T](d)). The two absorption mechanisms lead to different regimes 
of frequency redistribution of the reemitted photons. 
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C. Emitted photon 


As mentioned before, we consider a low-pressure vapor where atomic collisions are not frequent. In such 


a situation, and considering the low intensity of the incident radiation 29|, the scattering is elastic in the 


atomic rest frame, i.e. the emitted frequency is the same as the absorbed frequency 


li- 


In the laboratory frame, the detuning Si of the emitted photon is Doppler-shifted, Si = SA,i + where 
the direction of the emitted photon, of unit vector n', is drawn from an isotropic distribution. The Doppler 
shift {Si — Si-i) acquired in the ith emission constitutes the mechanism of frequency redistribution in the 
laboratory frame. If n is a unit vector in the direction of the incident photon, the atomic velocity component 
perpendicular to the plane (n, n') does not play any role in the frequency change is this scattering event and 
the emitted frequency, as measured in the laboratory frame, is given by: 

Si = Si-i + — (Vj| (cosd — 1) -I- V±sin0^ , (2) 

where Vj_ is the velocity component perpendicular to n, in the (n, n') plane and 9 is the angle between ft 
and n'. The Doppler frequency shift resulting from the isotropic reemission process therefore involves two 
of the atomic velocity components, Vj| and V±. 


D. Boundaries of the system 

We consider an atomic vapor with no boundaries (we call it infinite vapor). Simulating light transport in 
such an infinite vapor helps us understand the frequency redistribution mechanism at each scattering event. 
Furthermore, we can observe the probability distribution of the photons step length in order to get a deeper 
insight into the statistically anomalous properties of radiation trapping. 
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III. RESULTS 


A. Frequency redistribution 


1. First scattering event 


In this section we discuss the results of the simulation of radiation trapping in an infinite vapor. We plot 
in Figure [2] the PDF of the photons detuning after the first scattering event. The frequency is measured in 
the laboratory frame for several initial excitation frequencies (Sj). One may observe that for small excitation 
detunings, the reemission occurs around the line center while for excitation far from resonance the emission 
is centered at the incident frequency. Those curves are equivalent to the redistribution function R{di,6i) 


given by Hummer 
in Figure [H 


15l |. Curves calculated using equation 3.12.1 given Hummer 


15j| are shown as solid lines 


The qualitatively different behaviors of “into resonance” (|(5/| < S^) and “wing” excitation (jd/l > Sl) are 
clearly shown in Figure [5] They are directly related to the two distinct mechanisms of absorption discussed 
in section Hi HI For slightly-detuned excitation, the photon is absorbed by an atom with Vj| ss X6i, i.e., the 
photon frequency is Doppler-shifted to the line center in the atomic rest frame {6aj ~ 0). Reemission also 
occurs at <5^,1 ~ 0 and the frequency is Doppler-shifted in the laboratory frame. For excitation at <5/ = 0, 
Psi{V\\) is centered at zero velocity. The contribution of the very narrow distribution of the parallel compo¬ 
nent to the speed value is negligible and the PDF of the speed of the absorbing atoms approximately follows 
a 2D Maxwell-Boltzmann distribution (see Figure El). The most probable speed for the 2D MB distribution 
is smaller than that for the 3D MB distribution, resulting in a Doppler-shifted reemission narrower than 
the absorption profile [30| and centered at di = 0. For excitation detunings smaller than 5l {Sl = I.TFd 


in our system), the absorbing atom has a well dehned velocity component parallel to the incoming photon 
Vll Ri X6i and a distribution of velocity component normal to the incoming photon centered at zero. The 
atom speed therefore cannot be smaller than | V|| | (see Figure El and the acquired Doppler shift in emission 
is essentially given by (V||/A) (cos0 — 1) « Si{cos6 — 1) (see EqlJ). For an isotropic emission, cos{9) is 
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FIG. 2: (Color online) Probability density function for the frequency emitted after the first scattering event (<5i). The arrows 
in the Figure indicate the excitation frequency (Sj) for each curve. Calculated R{5i,5i) using equation 3.12.1 given by 

Hummer lal are shown as solid lines. 


uniformly distributed between —1 and 1, resulting in a plateau-like emission spectrum, of width ~ 2Si and 
amplitude proportional to 1/Si. 


For excitations with a detuning larger than Sl, absorption occurs off resonance in the frame of atoms 
following a thermal 3D MB distribution (see Figure [3]). Reemission occurs at the absorbed frequency in the 
atomic frame {Sa,i around Si). In the laboratory frame the emission frequency distribution is centered at Si 
and has a Doppler width. 
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Speed / u 


FIG. 3: (Color online) Probability density function of the speed of the absorbing atoms, tor different excitation detunings Sj. 
The magenta solid line superposed to the 5/ = 0 curve is a 2D MB distribution at 294 K. The green solid line superposed to 

the Sj = —2r£) curve is a 3D MB distribution at 294 K. 


2. Many scattering events 

We now turn to the observation of how the distribution of the emitted frequency evolves with the number 
of scattering events. In Figured is shown the PDF for an excitation at line center (<5/ = 0) and in Figure [S] 
for excitation at 5i = —F £>. For (5/^0 but smaller than 5L^ the first diffusion produces an almost uniform 
distribution of width 25 1 around the line center, as discussed in section PlI A II As the emitted frequency 
after the first scattering obeys |(5i| < 5l^ each ensemble of photons with a given will produce after the 
second scattering a plateau of height oc \/5i and width 25i. This results in a PDF centered at <52 = 0 with 
width larger than F o and smaller than 26 j. The following scattering events are very close to the absorption 
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FIG. 4: (Color online) Probability density function of the emitted frequency for different scattering event number i (denoted 
by the numbers in the figure). The initial excitation is at the line center (5/ = 0). 


profile of the vapor, converging to a CFR situation. 


Let us examine the emitted frequency PDF for large excitation detunings (5i > 5l)- The memory of 
the incident frequency is partially kept after the first scattering (see Figure IH]). However, the Doppler 
broadening of the emission frequency (5i around the excitation frequency 5i results in some photons close 
enough to resonance (Ji < 5l) to be subsequently absorbed by atoms with = A(5i. Those photons are 
thus reeinitted with detuning <52 around the line center. As a result, for this second scattering, a double-peak 
appears in the frequency distribution: one peak is centered at S 2 = Sj and the other one centered at S 2 = 0. 
In the following scattering events (i > 1), the peak around the line center gets higher until it dominates 
the emission frequency PDF. This behavior results in a net diffusion of the emitted frequency of the whole 
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fig. 5: (Color online) Probability density function of the emitted frequency for different scattering event number i (denoted 

by the numbers in the figure). The initial excitation is (5j = —F^j. 

sample and has consequences for the transport of photons in the vapor. 

To illustrate the frequency diffusion in the vapor we plot in Figure [7^ the emission frequency of a photon 
as a function of the number of scattering events for an excitation detuning 5i = —2r d for two different 
realizations of the MC simulation. One may see that at least a few scattering events are necessary for the 
emitted frequency to be inside the region |^i| < 5l where absorption preferentially happens for an atom with 
Vj| = X5i. Thereafter, emission occurs around the line center and the emitted frequency profile seems to 
converge to CFR. We plot in Figure [TJo the probability density function of the number of events necessary 
before emission occurs in the region |(5i| < 5l for the initial excitation detuning 5i = —2T£) used in Figure 
Eli. For this specific detuning, we see that approximately 55% of the photons have their frequency in the 
range |i5i| < 5l after the hrst scatttering. Flowever, a number of 15 scattering events is necessary to put 95% 
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FIG. 6: (Color online) Probability density function of the emitted frequency for different scattering event number i (denoted 
by the numbers in the figure). The initial excitation is at 5/ = —2r£). 


of the photons in this range. 


The spectral line shape of emission is important to interpret the behavior of the diffusion of photons in 
a resonant vapor. The frequency redistribution regimes usually considered are PFR after a single diffusion 


and CFR as the asymptotic result of many diffusion events 


24| . The question that naturally arises is how 


many cycles of absorption-reemission are necessary to obtain CRF. A way to analyze this question is to 
monitor the width of the emission spectrum as a function of the number of scattering events, as shown in 
Figure [5] for different detunings (d/ = 0, d/ = —F^i and 5i = —2r£)). We measure the full width at 1/e 
of its maximum in the emission spectra. For i5/ = 0 we see that the first scattering produces an emission 
frequency distribution that is narrower than the absorption spectrum in the vapor in accordance with the 
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FIG. 7: (Color online) a) Frequency diffusion during the random walk of a photon in a vapor as a function of the number of 
scattering events. Two photon realizations (red diamond and black triangle) are shown for 5j = —2r£). Dotted lines in (a) 
corresponds to :£Sl. b) PDF of the number of scattering events before |5i| < 5l for Sj = — 2r£) (left axis) and its cumulative 

probability (right axis). 


discussion in section UlI A II After the third scattering event one can consider that the width does not change 
much, being equal to r^i. This situation configures a CFR regime. Notice that for the Voigt parameter used 
(a = 0.01, Td r) the lineshape is dominated by Doppler profile [30| for small detunings and our numerical 


simulation sampling does not allow to differentiate the Doppler from the Voigt width. For Sj = —F/j the 
emission spectrum after the first scattering is very broad since it has the plateau-like shape discussed above. 
For the following scattering events (i > 2) the width gets narrower and reaches a stable value from the 
fourth scattering event on (i > 4). For Sj = —2Fd we have measured the width of the peak around the line 
center. For the first three scattering events this peak is not resolved and is not reported. We see that for 
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small numbers of events the frequency distribution is broad and gets narrower until its width stabilizes at 
r£), around the ninth scattering event. 



i 


FIG. 8: (Color online) Full width at 1/e of the maximum of the emitted frequency PDF as a function of the scattering event 
number for various excitation frequencies Sj. For Sj = —2r£i only the peak around the line center is reported. For the first 
three scattering events such peak is not resolved and is not reported here. 


B. Photon step-length distribution 


The wings of the emission and of the absorption spectral distributions play a crucial role in the transport 
of radiation through a vapor, since a photon emitted in the spectrum wings, i.e. far from resonance center, 
travels a longer path before being absorbed than if it is emitted close to the center of the resonance 9|. 
In particular, the asymptotic behavior of the step-length PDF {P{1)) dictates the statistics of the diffusion 
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23|,l2^. Measurements of P{1) in the CFR configuration require the preparation of a radiation with spectral 
profile identical to the absorption one, which can be achieved, e.g., through multiple scattering of laser 
radiation inside an atomic vapor prior to sending these photons to the measurement cell 2y, [27[. Notice, 
however, that to determine if CFR is achieved, the criterion of stability of the emitted spectrum width is 
not the most adequate since it gives little information on the spectrum wings. We show in Figures (|9]) to 
dm the evolution of P(/) as a function of the number of scattering events for the three values of Sj used 
in Figures 21 [5] and |6l For ^7 = 0 (Figure |9]) we see the stability of P{1) from the third scattering event 
on, indicating that three scattering events are enough to consider that the CFR regime is reached. For 
5i = —F 75 (Figure ITOl) . subtle changes of P{1) occur until the PDF stabilizes from the sixth scattering event 
on. For Sj = —2T 75 (Figure[TT]), at least nine scattering events are needed to achieve a stable P{1). 


For Sj = 0 (Figure 21), P{1) after the first scattering decays faster than expected for a diffusion of 

radiation with a Voigt profile into a Voigt absorption medium 2^ . This is due to the fact that the emission 
spectrum after the first scattering is narrower than the Voigt lineshape and thus closer to a monochromatic 
spectrum. For Sj = — F 77 (Figure fTUl). P{1) stabilizes with an asymptotic power law l~^ typical of a Doppler 


radiation incident in a va; 


Dor with Doppler absorption profile. This is expected for CFR in a vapor with Voigt 


absorption profile 23|, I 23 with Voigt coefficient a = 0.01 in the range of l/ao exhibited in Figure [TUI Our 
MC simulation has not reached the asymptotic range P{1) oc l~^-^ {l/cto > 10^) corresponding to very rare 
events far in the wings of the spectral distribution. For Sj = — 2 r 77 (Figure [IT]) . P{1) exhibits a power-like 
law for intermediate values of step length (see linear part of the curve in a log-log scale in Figure [TThl. For 
the asymptotic long steps, P{1) behaves rather as an exponential decay (see linear part of the curve in a 
mono-log scale in Figure [TTb l. The long steps are taken by photons emitted far from resonance in the broad 
peak around 6i = Sj = —2Tu, i.e., in the wings of the absorption profile. As the Voigt absorption profile 
falls as (5“^, the absorption coefficient around —2F 77 changes slowly in the far wings and P{1) approaches a 
behavior described by a Beer-Lambert law. This corresponds to a change in the statistical behavior of the 
light transport from superdifusive (for | 57 | < 6l) to normal diffusion (for |( 57 | > Sl) UM- The emission peak 


in the wing of the resonance around <5^ = ^7 = — 2 F 77 is broadened with the number of scattering events (see 
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FIG. 9: (Color online) Probability density function of the step length {P{1)) of the photon before being reabsorbed by the 
vapor for different scattering event number i (denoted as numbers in the figure). Excitation is at the line center Sj = 0. The 
step length is normalized by the absorption coefficient at the line center oq- 

Figure [6l). As a result, for large numbers of scattering events, an almost constant emission profile appears 
in the wing of the resonance line. The consequence of the subsistence of this constant profile is that CFR is 
never reached, that is, the emission profile does not completely converge to the Voigt absorption profile and 
P{1) maintains its exponential-like asymptotic regime over large numbers of scattering events. 


IV. CONCLUSION 


We have developed a Monte-Carlo simulation to analyze the diffusion of light propagating in a resonant 
atomic vapor. For the first scattering event, our simulations give results similar to the theoretical approach 
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FIG. 10: (Color online) Probability density function of the step length {P(l)) of the photon before being reabsorbed by the 
vapor for different scattering event number i (denoted as numbers in the figure). Excitation is at Sj = —Td. The step length 

is normalized by the absorption coefficient at the line center qq. 


of the redistribution function developed by Unno 


14j and Hummer 


15|. We have obtained the evolution 


of the PDF of the emitted frequency after multiple scattering events and analysed its convergence to a 
Complete Frequency Redistribution. Using only the width of the emitted spectrum as a criterion to check 
CFR, we observe that an average of three scattering events is necessary to reach CFR when the excitation 
is tuned to the line center. On the other hand, an average of six scattering events is necessary to fully 
redistribute an excitation detuned one Doppler width away from the resonance; Nine scattering events are 
necessary for excitation at two Doppler widths from resonance. A second criterion can be used to analyse 
CFR: the Probability Density Function of the step length of photons between an emission and the subsequent 
absorption in the vapor, P{1). For excitation close to resonance |(5/| < Sl, the long tail of P{1) behaves as 
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FIG. 11: (Color online) Probability density function of the step length {P(l)) of the photon before being reabsorbed by the 
vapor for different scattering event number i. The arrow indicate the sequence of curves with 2 = 1,2,4, 6, 8, 9,10. Excitation 
is at 5/ = — 2r£). The step length is normalized by the absorption coefficient at the line center qq. 


a power law that changes its coefficient at each scattering event. For those frequencies the photons undergo 
a superdiffusive transport inside the vapor. However, for excitation far from resonance (|(5/| > Sl) P{1) 
exhibits an exponential decay typical of normal diffusive transport. The exponential behavior originates from 
the photons emitted around the excitation frequency. Such a far-from-resonance emission subsists for a large 
number of scattering events preventing CFR to be reached and the transport maintains normal diffusion 
characteristics. We believe that these MC simulations allow to better interpret data from experiments of 
light diffusion in a resonant vapor. A natural prospect of this work is to include in the simulations the 
multilevel structure of the Rb isotopes. Indeed, optical pumping can occur between different ground states 
of an alkali atom and, moreover, transfer of radiation can occur between the two Rb isotopes due, e.q., to 
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the proximity of the levels ®®55'i/2(i^ = 3) and ^'^5Si/2{F = 2). 
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